Charge ordering and magnetism in quarter-filled Hubbard-Holstein model 
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We study a two-dimensional Hubbard-Holstein model with phonons treated in the adiabatic limit. 
A Hartree-Fock decomposition is employed for the Hubbard term. A range of electronic densities 
are discussed with special emphasis on the quarter-filling (n — 0.5). We argue that the quarter- 
filled system is relevant for the electronic properties observed at the interface between LaAlOs and 
SrTiOs, where half-electron per unit cell is transferred to the Ti02 layer as a consequence of the 
polar discontinuity at the interface. In addition to presenting the overall phase diagrams, we identify 
an interesting charge-ordered antiferromagnetic phase for n = 0.5, which was also reported recently 
in the ab-initio study of the LaAlOg-SrTiOs interface. 
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I. INTRODUCTION 

The interfaces between bulk insulating oxides such as 
LaAlOs (LAO) and SrTiOa (STO) have recently become 
a topic of very active research [J, [3, d, 0, Si- The in- 
terest in such interfaces was triggered by the oberva- 
tion of unusually high in-plane conductivity at the Ti02- 
terminated interface 0. Materials with such high mo- 
bilities can find numerous applications in various fields 
of industry and can also be utilized to make new devices 
d, Q ■ The origin of this effect is beheved to be an elec- 
tronic reconstruction caused by the polar discontinuity 
at the interface Q. The conceptual idea is that in or- 
der to avoid a divergence of electrostatic potential, 0.5 
electrons per unit cell are transferred to the Ti02 layer. 
These electrons then behave as a quasi two-dimensional 
electron gas leading to large mobilities. 

Subsequently, experiments were carried out at variable 
oxygen pressure and a strong dependence in the transport 
measurements was reported Q • In the presence of high 
oxygen pressure, an insulating behavior in the resistivity 
was observed at low temperatures. This indicates that 
perhaps oxygen vacancies play a crucial role in the exis- 
tence of high interfacial conductivity and the system may 
actually be insulating in the absence of oxygen vacancies. 
In addition, an external magnetic field has a strong ef- 
fect on transport 01- Applying magnetic field leads to a 
gain in the conductivity at low temperatures. An expla- 
nation for the insulating behavior of the resistivity was 
suggested to be connected to the magnetism and perhaps 
to the Kondo effect. 

While the concept of electronic charge reconstruction 
or oxygen vacancies can provide a reasonable explanation 
for the observed conductivity, there is no rigerous under- 
standing for the origin of magnetism at the interface. 
Hence, the understanding of magneto-transport is also 
rather incomplete. The present theoretical understand- 



ing of the interface properties has been largely based on 
the Density Functional Theory (DFT) calculations. Ex- 
istence of a charge ordered state has been pridicted by 
these calculations Q . At the level of a minimal model, 
a charge-ordered state can be obtained within the ex- 
tended Hubbard model [§]. But, such a state should 
be non-magnetic, suggesting that there might be an al- 
ternate mechanism active in these systems which leads 
to the charge ordering. Recent LDA-I-U studies of the 
LAO/STO interface have shown that structural distor- 
tions are present in the vicinity of the interface indicating 
the presence of an electron-lattice coupling. In this study 
a charge-ordered antiferromagnetic state was found to be 
the groundstate fld^ . 

In this paper we study the two-dimensional Hubbard- 
Holstein model as a simplest model capturing the effects 
of both, the electron-electron and the electron-lattice 
interactions. We find that the magnetic moments are 
formed as a consequence of polaron formation. These 
magnetic moments are found to be antiferromagneti- 
cally correlated at finite densities. At quarter filling, the 
groundstate is antiferromagnetic and charge-ordered in 
agreement with the findings of recent LDA-I-U calcula- 
tions. 



II. MODEL AND METHOD 

We consider a one-band Hubbard-Holstein model on a 
square lattice with the Hamiltonian: 
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Here, c.^^ and cj^ are the electron annihilation and cre- 
ation operators with electronic spin a ="1,1- n-i — 
TT-ii + riii is the charge density operator at site i with 

nia- = cl^c^a- The average electronic charge density is 
denoted by n and Xi denote the volume contraction and 
expansion of the oxygen octahedra, which couples to the 
variations in the charge density. The hopping parameter 
t is set to 1, therefore all energy scales arc in units of 
t. U is the strength of on-site Hubbard repulsion, and A 
is that of electron-lattice coupling. The lattice stiffness 
constant K is set to 1. 

In the present study the lattice distortions Xi arc 
treated in the adiabatic limit. In the absence of the Hub- 
bard term, the electronic Hamiltonian is bilinear in an- 
nihilation and creation operators, with the background 
potential provided by the configuration {xi\ of lattice 
distortions. The groundstate therefore, corresponds to 
the lattice configuration which minimizes the total en- 
ergy. In the absence of the kinetic energy term [t = 0), 
the problem reduces to N replicas of the one on a single 
site, where N is the number of lattice sites. The total 
energy for a single site is given by, 

E = -A.T,((n,) -n)+a;2/2. (2) 

Here and below {A) denotes expectation value of the op- 
erator A. Minimization of the energy E w.r.t the classical 
variable x^ leads to Xi = X{{ni) — n). For finite i, how- 
ever, the kinetic energy term also contributes to the total 
energy, 

-XJ2^^i{n^)-n) + l/2J2^^■ (3) 

i i 

Note that Xi is not a specified potential, but has to be 
determined self-consistently with the distribution of the 
electronic charge density. We compute the self-consistent 
lattice potential in the following scheme: Start with an 
arbitrary configuration of lattice distortions {xi}. Diag- 
onalize the Hamiltonian to generate the eigenvalues and 
eigenvectors. Compute the electronic charge density (rii). 
Use the relation Xi — X{{ni) — n) at each site to generate 
the new configuration for {xi}. Repeat the process un- 
til the old and new charge densities match within given 
error bar. 

In order to include the Hubbard term into the above 
self-consistent formalism, we treat the Hubbard term 
within Hartree-Fock approximation. The Hartree-Fock 
decomposition of the Hubbard term leads to, 

Hu = U^{n^'^)n^i + n,^{n,i) - (ri,-f) (ri,x). (4) 

i 

Now the self-consistency cycle requires the convergence of 
(nil) and (riii) individually. The generic problem with 



the self-consistent methods is that they need not lead 
to the minimum energy solution. Therefore, we use a 
variety of ordered and random initial states for the self- 
consistency loop and select the converged solution with 
the lowest energy. 

The Hubbard-Holstein model contains a variety of in- 
teresting phases and phenomena including, superconduc- 
tivity, charge- and spin-density wave formation, phase 
separation, and polaron and bipolaron formation. For 
this reason, the Hubbard-Holstein model has always been 
of interest in different contexts [ill, 113, El, [11, [ll, [11. 
Most of the earlier studies on this model were focused at 
or near half- filling. The quarter- filled case has not been 
analyzed in much detail except in one-dimension [l7j . 



III. RESULTS 
A. Dilute limit 

We begin by analyzing the case with very few electrons. 
Consider the Hamiltonian Eq. (1) with a single electron. 
The Hubbard term is inactive and for small A the ground- 
state wavcfunction corresponds to a Bloch wave. Since 
the lattice remains undistorted i.e. Xi = 0, the only con- 
tribution to the total energy is from the kinetic energy. 
For a single electron in a 2D square lattice, the lowest 
eigenvalue is — 4t, which is also equal to the total kinetic 
energy. Upon increasing the value of A, the energy is 
gained via the Holstein coupling term by self-trapping of 
the electron into a single polaron (SP). In mean field the 
trapping occures only when the energy of the SP state is 
lower than —At. Within a simple analysis, where we as- 
sume an ideal trapping of the electron at a single site, Eq. 
(2) leads to, Esp = -A^ + \^ /2 = /2. Hence, the 
critical value of A required for trapping a single electron 
into a single polaron is given by Af = 2\/2. 

Now consider the case of two electrons. In addition 
to the possibility of trapping the electrons as two single 
polarons, it is also possible to find a bipolaron (BP) so- 
lution. Infact the BP solution has a lower energy than 
two SPs and the critial value of A required to form a BP 
is given by X^^ — 2. The tendency to form bipolarons 
is clearly suppressed by the repulsive energy cost of the 
Hubbard term. From a very simple analysis of the two 
electron case one obtains: Efree ~ —8, Ebp ~ —2X^ + U 
and Esp ~ — A^. Looking for various energy crossings 
as a function of A and U one obtains the phase diagram 
shown in Fig. 1(a) for the three phases considered above. 
The solid lines arc from the simple analysis described 
above and the symbols represent the boundary values 
from the self-consistent numerical calculation. Note that 
this phase diagram corresponds to the case of two elec- 
trons in an infinite lattice and therefore refers to n — > 
in terms of fractional electronic filling in the thermody- 
namic limit. 

Assuming that U is large so that we are in the regime 
of single polaron formation, we estimate the effective 
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FIG. 1: (a) Phase diagram in the hmit of low electron den- 
sity in the parameter space of electron-lattice coupling A and 
the Hubbard repulsion U . Single polaronic and bipolaronic 
regimes are denoted by SP and BP, respectively. The solid 
lines are from a strong-trapping analysis (see text) and the 
symbols are the results of numerical calculations, (b) Total 
energy as a function of the distance between two single po- 
larons. The circles (squares) are for the parallel (antiparallel) 
spins of the single polarons. Inset shows the effective magnetic 
coupling between two self-trapped electrons as a function of 
the distance between them. 



interaction between two single polarons by calculating 
the total energy as a function of the distance between 
them. The energy difference between the spin-alligned 
and spin-antialligned single polarons provides an esti- 
mate for effective magnetic interaction between two po- 
larons. The energy variations are shown in Fig. 1(b), 
suggesting a repulsive and antiferromagnetic interaction 
between the localized magnetic moments. The energy 
difference AE' = E^^ — E^^ is plotted in inset in Fig. 
1(b). Positive values of Ai? for all R show that the two 
trapped moments prefer to be antiferromagnetic for all 
distances. In fact, the strength of the interaction is al- 
most vanishingly small for R > 2, suggesting the absence 
of any ordered magnetic state for low densities. We will 
see in the following sections that the above analysis of 
the dilute limit provides a very simple understanding of 
the phases that occur at finite densities, in addition to 
clarifying the basic competing tendencies present in the 
Hubbard-Holstein model. 



B. Generic Electron Densities 

For analyzing the system at higher electron densities 
we employ the self-consistent method described in the 
previous section. For the converged solution with mini- 
mum energy, we compute the charge structure factor, 

A.(q) = A^-2^((n,) -n)((n,) -n) e-i' 

ij 

and the spin structure factor. 
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FIG. 2: (Color online) The charge and spin structure factors 
at various q as a function of average charge density n aXU = 6 
for, (a) A = and (b) A = 2. The total density of states for 
selected n axe shown for, (c) A = and (d) A = 2. The dotted 
(blue) curves in (c) and (d) are shifted along the y-axis for 
clarity. 
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with Si = [riif — ni^)/2. Various ordered phases are in- 
ferred from the peaks in these structure factors. Fig. 
2(a) shows the spin structure factor at q = (tt, tt), which 
is a measure of antiferromagnetic correlations, and at 
q = (0, 0), which is indicative of a ferromagnetic behav- 
ior. At A = The system is antiferromagnetic (AFM) 
at and near n = 1, becomes ferromagnetic (FM) for 
0.7 < n < 0.9, and eventually becomes paramagnetic 
(PM). The antiferromagnetism at half- filling arises as a 
consequence of the nesting feature of the Fermi surface. 
The ferromagnetism at intermediate densities can be un- 
derstood within a Stoner picture which suggests that the 
repulsive cost coming from the Hubbard term can be re- 
duced by a relative shift of the spin-up and spin-down 
bands. At A = 2, the antiferromagnetic regime near 
n = 1 broadens (see Fig. 2(b)). The ferromagnetism is 
absent. Near n = 0.5 we find peaks in the charge struc- 
ture factor at (tt, tt), which indicates a charge-ordered 
(CO) state. Simultaneous peaks are found in the spin 
structure factors at (0, tt) and (tt, 0) pointing towards 
the existence of a nontrivial state with simultaneous ex- 
istence of charge- and spin-ordering. All the results pre- 
sented in this paper are for N = 32^, the stability of these 
results has been checked for system sizes up to iV = 40^. 
To further analyze the nature of electronic states we 
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FIG. 3: (Color online) The U-n phase diagrams for the 
Hubbard-Holstein model for (a) A = 0, and (b) A = 2. Ferro- 
magnetic (FM), antiferromagnetic (AF), paramagnetic (PM) 
phases are present in case of the pure Hubbard model (A = 0). 
Single polaron (SP), bipolaron (BP) and charge ordered an- 
tiferromagnetic phases also become stable for A = 2. 



compute the total density of states (DOS) as, 

j/tt 



N{lj) = N-^ J2 ^("^ - ~ J2 



72 + (w - e,;)2 



Here, a denote the eigenenergies corresponding to the 
minimum energy configuration. The (5-function is ap- 
proximated by a Lorentzian with width 7. We use 
7 = 0.04 in the calculations. A clean gap in the DOS 
is observed only for n=l in the absence of A (see Fig. 
2(c)). In the FM regime, a two-peak structure represents 
a shifted spin-up and spin-down band, which is consis- 
tent with the Stoner picture of magnetism in Hubbard 
model. Eventually at low-density the DOS begins to re- 
semble the free electron tight-binding DOS. More inter- 
esting features are observed in the DOS at A = 2 shown 
in fig. 2(d). The clean gap originating from the AFM 
state survives down to n '--^ 0.85. The gap opens up once 
again at quarter-filling (n = 0.5). This correlates per- 
fectly with the signatures found in the structure factor 
calculations shown in Fig. 2(b). 

The results for various t/ at A = and A = 2 are 
summarized into two phase diagrams. The U-n phase di- 
agram for A = is shown in Fig 3(a). Antiferromagnetic, 
ferromagnetic and paramagnetic states are found to be 
stable in agreement with prev ious results on the Hubbard 
model in two-dimensions ]l8l.[l9l.[20t. Fig. 3(b) shows the 
U-n phase diagram for A = 2. For low U, the system be- 
comes a bipolaronic insulator. Although the Xf^ ^ 2 
for a single BP, it can be much lower for finite density, 
suggesting that it is easier to trap many bipolarons as 
compared to a single BP [2l[. The charge ordered state 
at half-filling can be simply viewed as a checkerboard 
arrangement of bipolarons, although the concept of an 
isolated BP does not really hold for such large densities. 
The charge-ordered state exists even below the critical 
coupling required for BP formation due to the nesting 
feature present in the Fermi surface at half filling. The 



FIG. 4: A — f/ phase diagrams at (a) half filling, and (b) 
quarter filling. 



half-filled CO state undergoes a transition to an antifer- 
romagnetic state near [/ = 4. The region of antiferro- 
magnetism grows with increasing U in contrast to the 
pure Hubbard model. The PM state still exists for small 
A but the FM state is absent. A large region of phase 
space is taken by the single polaronic state for large U. 
No magnetism is found at low densities, since these sin- 
gle polarons are magnetically non-interacting due to the 
large inter-polaronic separations. At large densities, how- 
ever, There are antiferromagnetic correlations between 
these single polarons. This is consistent with the effective 
magnetic interactions found between two single polarons 
(see Fig. 1(b)). These effective antiferromagnetic inter- 
actions are the origin of the growth in the AFM regime 
near n = 1. 



C. Half- and quarter-filling 

The Hubbard-Holstein model at and near half-filling 
has been studied previously [13, [H, . The existence 
of spin- and charge-density waves was reported. The 
possbility for an intermediate metallic phase was also re- 
ported in a one-dimensional model with dynamical effects 
for lattice [l^. Fig. 4(a) shows a,U — X phase diagram at 
half-filling. The system is either charge ordered or anti- 
ferromagnetic and therefore, the DOS is always gapped. 
The boundary separating the CO and the AFM states 
fits very well as \J=X^, which happens to be the bound- 
ary separating the SP and BP regime in the low density 
limit (see Fig. 1(a)). This suggests that the CO phase 
can be viewed as a checkerboard pattern of bipolarons, 
at least for large values of A. The origin of the CO or the 
AFM phase at small values of U and A is related to the 
existence of nesting in the Fermi surface with a nesting 
wavevector q = (tt, tt). 

The most interesting result of this paper is the obser- 
vation of a CO-AFM state at quarter filling. We plot 
the U — X phase diagram for quarter-filled system in Fig. 
4(b). Unlike the half-filled case, the small U, small A 
regime corresponds to free electron behavior. The SP 
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FIG. 5: (Color online) Real space patterns for charge density 
(upper row) and spin density (lower row) at i7 = 6 and A = 1 
(left column), and A = 2 (right column). Note that the spin 
state for A = 2 is a G-type antiferromagnet if one rotates the 
lattice by 45° and consider the square lattice of "occupied" 
sites. 



state is found to exist for large U and there is a large win- 
dow where a charge-ordered AFM state exists. We find 
that a self-consistent solution corresponding to a CO-FM 
state can exist only for U > 10, but it is still higher in 
energy than the CO- AFM. For U < 10, the CO-FM state 
is not stable and therefore, the charge ordering occures 
only when it is accompanied by an AFM ordering. This 
leads to a very interesting implication for the effects of 
external magnetic field. Enforcing ferromagnetism by ap- 
plying an external magnetic field to the CO-AFM state 
would lead to a melting of the charge order and hence to 
a large negative magnetoresistence. 

To further investigate this unusual state at quarter- 
filling, we show the real-space data for charge density rii 
and spin density Si in Fig. 5. A weak charge-ordering 
is already present at A = 1 (see upper-left panel). This 
charge ordering is accompanied by a stripe-like spin or- 
dering, where spins are arranged ferromagnetically along 
one direction and antiferromagnetically along the other 
(lower-left panel). A very clear CO pattern emerges 
for the larger value of A (upper-right panel), which oc- 
cures together with an AFM arrangement of the spins 
(lower-right panel). In a strong coupling scenario, it is 
easy to understand the CO-AFM state within the picture 
of effective magnetic interaction between single polarons 



presented in Fig. 1(b). Assuming that the "occupied" 
sites in the charge ordered state can be viewed as single 
polarons, an effective antiferromagnetic interaction be- 
tween them is strongest at distance ^/2, hence leading to 
a magnetic structure which is a G-type AFM order for 
the square lattice constructed out of the " occupied" sites 
only. For smaller values of A the charge disproportion- 
ation in the CO state is much smaller and the effective 
magnetic interaction picture can not be pushed to this 
weak coupling regime. 



IV. CONCLUSIONS 

We have presented groundstate properties of the 
Hubbard-Holstein model in two dimensions in the adi- 
abatic limit for the lattice distortions. We use a self con- 
sistent method for generating the static lattice configu- 
rations in combination with a Hartree-Fock decoupling 
of the Hubbard term. Interestingly, the charge-ordered 
antiferromagnetic state that we find at quarter filling was 
shown to be the groundstate for the LAO/STO interface 
in recent DFT calculations [l3|. Within our analysis the 
charge ordering in this state occures only in combina- 
tion with the AFM ordering, as we find that the CO-FM 
state is unstable. Therefore, the charge ordering could 
be melted by applying an external magnetic field, lead- 
ing to a large negative magnetoresistence. We argue that 
this model is relevant for the LAO/STO interface since, 
(i) it provides a possibility for the formation of magnetic 
moments, (ii) leads to a CO-AFM groundstate in agree- 
ment with the recent LDA+U studies, and (iii) contains 
the possibility for large negative magnetoresistence via 
a magnetic-field-induced melting of the charge ordered 
state. 
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